Superconducting transition temperature in thin films 
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By considering the Ginzburg-Landau model, compactified in one of the spatial dimensions, and 
using a modified Matsubara formalism, we determine the dependence of the superconducting tran- 
sition temperature (T c ) of a film as a function of its thickness (L). We show that T c is a decreasing 
linear function of as has been found experimentally. The critical thickness for the suppression 
of superconductivity is expressed in terms of the Ginzburg-Landau parameters. 
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In last decades, a large amount of work has been done 
on the Ginzburg-Landau model applied to the study of 
the superconducting transition, both in the single com- 
ponent and in the N- component versions of the model, 
using the renormalization group approach. The state of 
the subject, for type-I and type-II superconductors and 
related topics, can be found for instance in Refs. ]l]-f§. 
In another related topic of investigation, there are sys- 
tems that present domain walls as defects, created for in- 
stance in the process of crystal growth by some prepared 
circumstances. At the level of effective field theories, in 
many cases, this can be modeled by considering a Dirac 
fermionic field whose mass changes sign as it crosses the 
defect, meaning that the domain wall plays the role of 
a critical boundary separating two different states of the 
system |7],|8|]. Questions concerning stability and the ex- 
istence of phase transitions may also be raised if one con- 
siders the behavior of field theories as function of spacial 
boundaries. Studies on confined field theory have been 
done in the literature since a long time ago. In particu- 
lar, an analysis of the renormalization group in finite size 
geometries can be found in Refs. ||Jl(J. These studies 
have been performed to take into account boundary ef- 
fects on scaling laws. The existence of phase transitions, 
would be in this case associated to some spatial param- 
eters describing the breaking of translational invariance, 
for instance the distance L between planes confining the 
system. In this situation, for Euclidean field theories 
the Matsubara formalism applies for the breaking of in- 
variance along any one of the spacial directions. Studies 
of this type have been recently performed fll]Jl^] , con- 
cerning with the spontaneous symmetry breaking in the 
A</) 4 theory. In particular, if one considers the Ginzburg- 
Landau model confined between two parallel planes, thus 
describing a superconducting film, the question of how 
the critical temperature depends on the thickness L of 
the film can be raised. 

Under the assumption that information about general 
features of the behavior of superconductors, in absence 
of magnetic fields, can be obtained in the approximation 
which neglects gauge field contributions in the Ginzburg- 



Landau model, in this letter we examine this model with 
an approach different from the renormalization group 
analysis. We consider the system confined between two 
parallel planes and we use the formalism developed in 
Refs. p|, 12| to investigate how the critical temperature 
is affected by the presence of boundaries. From a physical 
point of view, we investigate how the critical temperature 
of a superconducting film depends on its thickness. 

We start with the Ginzburg-Landau Hamiltonian den- 
sity in the Euclidean D -dimensional space, in absence of 
magnetic fields, given by (in units with K = 1) 
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where A is the (renormalizcd) self-coupling constant, with 
the "bare mass" given by m\ = a(T — To), To being 
the bulk transition temperature of the superconductor 
and a > 0. We consider the system confined between 
two parallel planes, normal to the distance L 

apart from one another and use Cartesian coordinates 
r = (,x, z), where z is a (D — l)-dimcnsional vector, with 
corresponding momenta k = (fc,j,,q), q being a (D — 1) 
-dimensional vector in momenta space. The partition 
function is written as, 

Z = J XV*XVexp f-jf dx J dp- 1 * H{\ip\, |V<^ , 

(2) 

with the field ip(x,z) satisfying the condition of confine- 
ment along the x-axis, ip(x < 0, z) = ip(x > L, z) = 0. 
Then the field should have a mixed series-integral Fourier 
representation of the form, 
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where oj n = 2itn/L and the coefficients c„ and 6(q) cor- 
respond respectively to the Fourier series representation 
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over x and to the Fourier integral representation over the 
(D — l)-dimensional z -space. The above conditions of 
confinement of the x-dependence of the field to a seg- 
ment of length L allow us to proceed, with respect to 
the x -coordinate, in a manner analogous as it is done in 
the imaginary-time Matsubara formalism in field theory 
and, accordingly, the Feynman rules should be modified 
following the prescription 
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We emphasize that we are considering an Euclidean field 
theory in D purely spatial dimensions, so we are not work- 
ing in the framework of finite temperature field theory. 
Here, the temperature is introduced in the mass term of 
the Hamiltonian by means of the usual Ginzburg-Landau 
recipe. 

To continue, we use some one-loop results described 
in pj]Jl^| , adapted to our present situation. These re- 
sults have been obtained by the concurrent use of di- 
mensional and zeta-function analytic regularizations, to 
evaluate formally the integral over the continuous mo- 
menta and the summation over the Matsubara frequen- 
cies. We get sums of polar (L -independent) terms plus 
L-dependent analytic corrections. Renormalized quanti- 
ties are obtained by subtraction of the divergent (polar) 
terms appearing in the quantities obtained by applica- 
tion of the modified Feynman rules (Matsubara prescrip- 
tion) and dimensional regularization formulas. These po- 
lar terms are proportional to r-functions having the di- 
mension D in the argument and correspond to the in- 
troduction of counter-terms in the original Hamiltonian 
density. In order to have a coherent procedure in any 
dimension, these subtractions should be performed even 
for those values of the dimension D where no poles of T 
-functions are present. In these cases a finite renormal- 
ization is performed. 

In the following, to deal with dimensionless quantities 
in the regularization procedure, we introduce parameters 
c 2 = ™ 2 /47T 2 /i 2 , a = (L[i)~ 2 , g = 3A/87T 2 and 4> = tpo/fJ,, 
where ipo is the normalized vacuum expectation value of 
the field (the classical field) and fi is a mass scale. In 
terms of these parameters, the one-loop contribution to 
effective potential, adapted to the situation under study, 
is given by the well known expression Mm 
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where m (entering in c) is the renormalized mass for 
L = 00. Performing the Matsubara replacement (Q), 
the boundary-dependent (L-dependent) one-loop contri- 
bution to the effective potential can be written in the 
form 
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Now, using the dimensional regularization formula, 
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Eq. (H) reduces to 
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where f(D, s) is a function proportional to T(s — 

and zf (s — ; a) is one of the Epstein-Hurwitz zeta 
-functions, defined by 
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valid for Re{u) > K/2 (in our case Re(s) > D/2). 

The Epstein-Hurwitz zeia-function can be extended to 
the whole complex s -plane and we obtain, after some ma- 
nipulations p5[, the one- loop correction to the effective 
potential, 
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and K v are the Bessel functions of the third kind. 

Note that since we are using dimensional regulariza- 
tion techniques, there is implicit in the above formulas 
a factor fi 4 ~ D in the definition of the coupling constant. 
In what follows we make explicit this factor, the sym- 
bol A standing for the dimensionless coupling parameter 
(which coincides with the physical coupling constant in 
D — 4). We work in the approximation of neglecting the 
L-dependence of the coupling constant, that is we take A 
as the renormalized coupling constant. In this case, it is 
enough for us to use only one renormalization condition, 
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Since we are using a modified minimal subtraction 
scheme, where the mass (and coupling constant, if it is 
the case) counter-terms are poles at the physical values 
of s, the L-dependent correction to the mass is propor- 
tional to the regular part of the analytical extension of 
the Epstein-Hurwitz zeia-function in the neighborhood 
of the pole at s = 1. Thus the L-dependent renormalized 
mass, at one-loop approximation, is given by 
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On the other hand, if we start in the ordered phase, 
the model exhibits spontaneous symmetry breaking, but 
for sufficiently small values of T _1 and L the symmetry 
is restored. We can define the critical curve C(T C ,L C ) 
as the curve in the T x L plane for which the inverse 
squared correlation length , £~ 2 (L, L, ipo), vanishes in the 
L-dependent gap equation j9), 
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where ipo is the normalized vacuum expectation value of 
the field (different from zero in the ordered phase). In 
the disordered phase, in particular in the neighborhood 
of the critical curve, ifQ vanishes and the gap equation 
reduces to a L-dependent Dyson- Schwinger equation, 
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After steps analogous to those leading from Eq.([|) to 
Eq.(O), Eq.((l3) can be written in the form 
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If we limit ourselves to the neighborhood of criticality, 
m 2 (T,L) ps 0, we may investigate the behavior of the 
system by using in Eq.(|l6|) an asymptotic formula for 
small values of the argument of Bessel functions, 



(*~0), 



which allows after some straightforward manipulations, 
to write Eq.(16) in the form 
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where ((D — 2) is the Riemann zeto-function, ((D — 2) = 
Y,n=i( l l nD ~ 2 )i defined for D> 3. Taking m 2 (T,L) = 
and m 2 , = a(T — To) m Eq.jT^), we obtain the critical 
curve in the T x L plane for Euclidean space dimension 
D (D> 3), 

a[Tc ~ + " ^ L c~°^ D 2 ) = °' ( 19 ) 

For D — 3 the Riemann zeia-function in Eq.(|l9|) has 
a pole. We can not obtain a critical curve in dimension 
D < 3 by a limiting procedure from Eq . (|l9|) . For D = 3, 
which corresponds to the physically interesting situation 
of the system confined between two parallel planes em- 
bedded in a 3-dimensional Euclidean space, Eq.fll9|) be- 
comes meaningless. To obtain a critical curve in D < 3, 
we perform an analytic continuation of the zeta-function 
to values of the argument z < 1, by means of the 
reflection property of zeta-functions 
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which defines a meromorphic function having only one 
simple pole at z = 1. For D = 3, a mass renormaliza- 
tion procedure can be done as follows: remembering the 
formula 
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where 7 = 0.57216 is the Euler constant, we define the 
renormalized mass 777 as 
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Taking this renormalized mass equal to zero leads to the 
critical curve in dimension D = 3, given by 
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In Eq.(23), Tq corresponds to the transition tempera- 
ture for the material in absence of boundaries (L c — > 00), 
that is, to the bulk transition temperature. We see then 
that, in a film made of the same material, the critical tem- 
perature is diminished by a quantity proportional to the 
inverse of its thickness. Also, we see that there is a min- 

(17) 

imal film thickness below which superconductive is 
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suppressed, which is given by (identifying the Ginzburg- 
Landau parameter (3 = A/i) 
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Such a linear dependence of T c with the inverse of the 
film thickness has been found experimentally in materials 
containing transition metals, for example, in Nb |17 1^] 
and in W-Re alloys ]20|] ; for these cases, it has been ex- 
plained in terms of proximity, localization and Coulomb- 
interaction effects. Notice that our result does not de- 
pend on microscopic details of the material involved nor 
accounts for the influence of manufacturing aspects, like 
the kind of substrate on which the film is deposited. In 
other words, the linear decreasing of T c as the film thick- 
ness is diminished emerges solely as a topological effect 
of the compactification of the Ginzburg-Landau model 
in one direction; other aspects, which may influence the 
transition temperature of the film, will show up experi- 
mentally as deviations of this linear behavior. 

Here, in a field theoretical framework, we have shown 
that quantum corrections to the mass in the Ginzburg- 
Landau model compactified in one of the spatial dimen- 
sions, in one loop-order, leads to the linear dependence of 
T c with the inverse of the thickness for a film, supercon- 
ductivity being suppressed at a minimum critical thick- 
ness One expects, however, that the inclusion of the 
L-dependence of the coupling constant in first-order may 
lead to a small correction to the linear behavior obtained. 
Our treatment can also be extended to consider external 
magnetic fields, but these issues will be discussed else- 
where. 
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